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The competition between Kondo screening and indirect magnetic exchange is studied for a system 
with geometrical frustration using dynamical mean-field theory (DMFT). We systematically scan 
the weak- to strong-coupling regime of the periodic Anderson model on the triangular lattice for 
a wide range of fillings n. The magnetic phase diagram is derived using a site-dependent DMFT 
approach by self-consistent mapping onto three independent single-impurity models corresponding 
to the three correlated / orbitals in the unit cell. At half-filling, the system is a non-magnetic 
Kondo insulator for all considered interaction strengths U > 0 which immediately develops into 
a non-magnetic metallic Kondo-singlet phase for fillings slightly below half-filling. On the other 
hand, indirect magnetic exchange between the / moments results in antiferromagnetic order at lower 
fillings. The antiferromagnetic and the Kondo-singlet phases are separated in the U-n phase diagram 
by an extended region of partial Kondo screening, i.e., a phase where the magnetic moment at one 
site in the unit cell is Kondo screened while the remaining two are coupled antiferromagnetically. 

At even lower fillings the system crosses over from a local-moment to a mixed-valence regime where 
the minimization of the kinetic energy in a strongly correlated system gives rise to a metallic and 
partially polarized ferromagnetic state. 

PACS numbers: 71.27.+a,75.10.Jm,75.30.Mb 


I. INTRODUCTION 

The periodic Anderson model is the model of choic e to 
describe heavy-fermion m ater ials realized in crystal^^^ 
or in quantum simulations.!^*^ It generically describes a 
band of light conduction electrons of bandwidth W hy¬ 
bridizing, with matrix element V, with a narrow band 
of / electrons located at energy e/. Since the / band 
is narrow, Coulomb correlations are important and are 
taken into account by an on-site Hubbard interaction U. 

The emergent many-body scales depend very much 
on the choice of bare parameters. In the absence of 
electronic correlations, hybridization leads to an intrin¬ 
sic /-band width F = ttV'^/W which gives an energy 
scale to assess the impact of the Hubbard U. In par¬ 
ticular, the formation of a local magnetic moment in 
a metallic host has been studied in Ref. HI It corre¬ 
sponds to a choice of bare parameters where the lower 
(upper) Hubbard band is below (above) the Fermi en¬ 
ergy /i, namely ef+r/2<fi<ef — r/2 + U. For strong 
U in the local-moment regime the effective low-energy 
physics can b e ap proximated by the more simple Kondo- 
lattice modelP*^ Here, charge fluctuations on the / or¬ 
bitals are completely suppressed, and a super-exchange¬ 
like mechanisrrP yields a magnetic energy scale, namely 
a local antiferromagnetic exchange J = 8V'^/U between 
the local and the conduction-electron moments. The lo¬ 
cal magnetic moment corresponds to a local Kramers 
doublet and is thereby - in the absence of correlations 
- protected by time-reversal symmetry. 

The residual entropy can be quenched by different 
and competing mechanisms. Magnetic ordering breaks 
time-reversal symmetry and is driven by the Ruderman- 


Kittel-Kasuya-Yosida (RKKY)P interaction. The cor¬ 
responding energy scale is set by the effective coupling 
strength Jrkky)?) = = 0) where is the 

conduction-electron spin susceptibility. It is an indirect 
interaction which is mediated via magnetic polarization 
of the conduction elec tron s. This energy scale competes 
with the Kondo scal^^ES! given by Tk oc . For 

temperatures below Tk the local magnetic moment is 
screened through the formation of a many-body entan¬ 
gled spin-singlet state with the conduction-electron spin 
degrees of freedom. The competition between RKKY 
coupling and Kondo screening leads to the famous Do- 
niach diagranP^* and to corresponding quantum phase 
transitions.fi^ 

Some heavy-fermion materials, such as CePdAl, are 
synthesized on frustrated geometries.*^^ This introduces 
another energy scale in the problem associated with the 
release of frustration via a mechanism of partial Kondo 
screening (PKS).^Here, a site-selective Kondo effect al¬ 
leviates the frustration thus allowing the remnant spins 
to order magnetically via the RKKY interaction. Such 
site-dependent screening can occur spontaneously or can 
reflect chemically different environments in compounds 
with large unit cells. 

The mechanism of partial Kondo screeni ng ha s at¬ 
tracted considerable attention in the past.E^E^ The 
purpose of the present paper is to study the effect 
in the Anders on la ttice beyond the static mean-field 
approximation.lA^liSI This is achieved by applying a vari¬ 
ant of the dynamical mean-field theory (DMFT))^ for 
the Anderson model on the triangular lattice where the 
different correlated orbitals in the unit cell are treated in¬ 
dependently, similar to a real-space DMFT approach.!^ 
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We will show that the competition between RKKY cou¬ 
pling and Kondo screening, supplemented by lattice frus¬ 
tration leads to a remarkably rich phase diagram includ¬ 
ing a PKS phase emerging in the local-moment regime 
at the border between paramagnetic heavy-fermion and 
magnetically ordered phases. 

The Anderson lattice has richer physics than the 
Kondo-lattice model since it allows for charge fluctua¬ 
tions on the / sites. In conjunction with strong spin-orbit 
coupling, for example, this natur ally le ads to the concept 
of a topological Kondo insu latoj^IJ^ which has argued 
to be realized in Here, we extend our study 

beyond the local-moment regime and consider magnetic 
phases in the mixed-valence regime. This is realized when 
the lower (or upper) Hubbard band overlaps with the 
chemical potential, i.e., if £/ -|- r/2 ~ /r. In this regime 
our numerical calculations indicate a completely different 
physics and predict ferromagnetic order in particular. We 
argue that the latter is reminiscent of itinerant-electron 
ferromagnetism caused by strong electron correlations in 
single-band models on frustrated geometries.^ 

Our study should be understood as a first step only 
which contributes towards a deeper understanding of the 
competition between magnetic order and Kondo screen¬ 
ing on frustrated lattice geometries: While the dynamical 
mean-field theory treats the local, temporal correlations 
exactly and correctly accounts for the Kondo effect, it 
also suffers from the simple mean-field-type description 
of spatial correlations. As concerns the effective RKKY 
interaction, the DMFT does capture its full spatial struc¬ 
ture and therewith the corresponding tendencies towards 
magnetic ordering but the feedback of non-local magnetic 
correlations on the one-particle Green’s function is ne¬ 
glected. We expect that this missing feedback will result 
in a somewhat biased description which probably overes¬ 
timates the instabilities against magnetic ordering. De¬ 
spite this and other possible deficiencies that are char¬ 
acteristic to any mean-field approach, we believe that 
the non-perturbative and internally consistent physical 
picture that is provided by the DMFT will serve as an 
important starting point for future studies, such as clus¬ 
ter and other extensions of the DMFT concept,!^ where 
the effects of short- and long-range correlations are pro¬ 
gressively included. Even this route cannot be expected 
to provide a final answer, given the complexity of the 
problem posed by strong correlations in fermionic models 
on two-dimensional frustrated lattices, and must be sup¬ 
plemented by complementary approaches, such as varia¬ 
tional wave functions (see Ref. [57] for an example). How¬ 
ever, the phase diagram derived from the site-dependent 
DMFT approach presented here will in any case provide 
a useful starting point and a valuable point of orientation 
in this general context. 

The article is organized as follows. The next section in¬ 
troduces the model, the site-dependent DMFT approach 
and the solver employed here. Section |HT] presents the 
numerical results. We discuss the DMFT phase diagram 
and analyze the different phases and mechanisms in de¬ 


tail. Conclusions are given in section IV 


II. MODEL AND METHOD 


The periodic Anderson model describes correlated “f” 
orbitals with a repulsive on-site interaction which locally 
hybridize with the “c” orbitals of a non-interacting sys¬ 
tem of itinerant conduction electrons. We study the 
Anderson model on the two-dimensional triangular lat¬ 
tice and consider a partitioning of the lattice into non¬ 
primitive unit cells containing three sites each, as shown 
in Fig. [2 Within the variant of the standard DMFT 
approach employed here, these sites will be treated as 
inequivalent (see below). Using standard notations, the 
Hamiltonian reads: 


H = 


H (r - r')c. 

rr' OiOc' cr 

+ V E 

roi(7 


U 


^ ^ fracrfraafra — afra — a ■ 


( 1 ) 


Here, r runs over the position vectors to the different 
unit cells, a G {A,B,C} refers to the sites within a unit 
cell, and a G {t,4-} is the spin projection, if^aa) 

creates an electron in the c (/) orbital with quantum 
numbers r,a,a. Conduction electrons are assumed to 
hop between nearest-neighboring sites, i.e., the hopping 
amplitude taa'i^ — r') = t ^ 0 ii r,a and r', a' are near¬ 
est neighbors. Furthermore, V is the local hybridization 
strength, and U is the strength of the Hubbard-type local 
interaction on the / orbitals. The one-particle energy of 
the / orbitals is £/ and for the c orbitals £c = tctQ(O). 

Dynamical mean-field theory (DMFT)P^ assumes 
that the self-energy on the / orbitals be local, 
UrQcr,r'a'cr'(w) = 6r,r'Saa''^<Ta-', and maps the lat¬ 
tice problem onto an effective single-impurity Ander¬ 
son model with one-particle parameters or, equivalently, 
with a hybridization function (w) that is deter¬ 
mined from the local element of the lattice Green’s func¬ 
tion Gioc,ct(t'(w) via the DMFT self-consistency condi¬ 
tion. This implicitly assumes that the dynamical mean- 
field Ao-cr'(w) is homogeneous. Consequently, only ho¬ 
mogeneous phases of the DMFT equations can be found 
in this way. In the real-space DMFT approaclP^ the 
self-energy is still assumed as completely local but inho¬ 
mogeneous solutions of arbitrary complexity are allowed 
by keeping the full spatial dependence of the local self¬ 
energy. fj/(^LU^ ,r'^act'^va’,aa'' Here, we 

employ a “site-dependent DMFT” by assuming that the 
local self-energy has possibly different elements on the 
different sites in a unit cell that is larger than a primi¬ 
tive cell. Otherwise, the self-energy is taken as homoge¬ 
neous: Erao-.T-'aV'(w) = 6r,r'3aa'Scra-'^aa-(^)- Restrict- 
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FIG. 1: (Color online) Periodic Anderson model on the tri¬ 
angular lattice. Right: Primitive unit cell (light gray, dashed 
lines) and unit cell (gray, solid lines), spanned by the vec¬ 
tors ai and a 2 , considered here. The latter contains three 
sites (A, B, C) treated independently within site-dependent 
dynamical mean-field theory. Left: For each site, a correlated 
(/) orbital with local interaction U couples to an uncorre¬ 
lated conduction-electron (c) orbital via the hybridization of 
strength V. The nearest-neighbor hopping t — 1 between 
conduction-electron orbitals sets the energy scale. 


ing ourselves to collinear magnetic phases for simplicity, 
we consider a possibly spin-dependent but spin-diagonal 
self-energy. 

For the above-mentioned partitioning of the triangular 
lattice, this ansatz for the self-energy means that the pe¬ 
riodic Anderson model is self-consistently mapped onto 
three independent and impurity models with possibly 
spin-dependent but spin-diagonal one-particle parame¬ 
ters. The impurity models can be solved independently 
but are coupled indirectly via the DMFT self-consistency 
equation. In particular, we do not impose any further 
condition on the spatial or spin dependence of Sq,o-(w). 
Thereby, we can account for different phases, in particu¬ 
lar for collinear magnetic phases, characterized by inho¬ 
mogeneous order parameters within a unit cell. 

Below, we briefly list the main equations of the site- 
dependent DMFT approach: Exploiting the remaining 
translational symmetry, Fourier transformation of the 
one-particle term of the Hamiltonian Q provides us with 
a 6 X 6 hopping matrix 


/£/ 

0 

0 

V 

0 


0 


0 

0 

V 

0 

0 

0 

£/ 

0 

0 

V 

V* 

0 

0 

£c 

£AB(k) 

eAc(fe) 

0 

V* 

0 

£ab(^) 

£c 

£Bc(fe) 

V 0 

0 

V* 

£ac(^) 

£bg(^) 

£c / 


for each wave vector k in the reduced Brillouin zone. We 
have 


eAB(fc) = t 

1 -1- 2 cos {ky/2) 

£Ac{k) = t 

1 _|_ g-*#fcx g-jjfey _(_ g-* 

eBc(fc) = t 

1 -1- -1- 


( 3 ) 


With this, and with a guess for the local but site- 
dependent / self-energy Sqct(w) (for a S {A,B,C}) we 
can start the DMFT self-consistency cycle by calculating 
the elements of the local lattice Green’s function via 


Gioc,7(5,ct (^) — ^ ^ ^ 


keBZ 


• + ti- e(k) - S^(w) 


, ( 4 ) 


7(5 


where 7 , <5 run over the 6 orbitals in the unit cell and 
where the 6 x 6 -matrix is diagonal and non¬ 

zero on the / orbitals only, /i is the chemical potential 
that is used to fix the total particle density. The local 
Green’s function is used to determine the hybridization 
functions of the three single-impurity Anderson models 
{a S {A, B, G}) as 


^acr{i^) = ^ + R- £f -^aa{0j) - — - . (5) 

^loc,aQ;,(7 (.^ ) 

Having defined the impurity models, the self-consistency 
cycle is closed by calculating the self-energy for 

each impurity model independently. 

The computational bottleneck of the DMFT cycle con¬ 
sists in the solution of the effective impurity problems. 
Here, we use the continuous-time quantum Monte-Carlo 
methocpS 122 l based on the hybridization expansion of the 
action of the respective impurity modeP^ at finite but 
low temperatures T. Since the interaction term is a 
density-density Hubbard-type interaction only, it i s ad- 
vantageous to employ the segment-picture variant .®EI] 
Following Ref. |311 this allows us to directly measure the 
impurity self-energy 'Eaa-i'ii^n) on the fermionic Matsub- 
ara frequencies iuJn- 


III. RESULTS 

DMFT calculations have been performed for the model 
Eq. Q with different chemical potentials /i to scan the 
interesting regime at and off half-filling n = 1 where n is 
given by 


a—A,B,C 

withni4^ = {flacrfrao) and Ua] = {cl^^Crao)- The Hub¬ 
bard interaction is scanned in the weak- to intermediate¬ 
coupling range 0 < C/ < 4 where the nearest-neighbor 
hopping t = 1 fixes the energy scale throughout the 
paper. Note that choosing t > 0 is convenient as this 
implies that the center of gravity of the total density 
of states (see Fig. is located close to the lower band 
edge. Symmetry-broken magnetic phases are therefore 
expected to occur for fillings below half-filling. We fur¬ 
thermore fix the hybridization strength at U = 1 and 
choose Ef = —U12 for the on-site energy of the / or¬ 
bitals. For strong U, this ensures that the occupancy 
of the / orbital at any site a in the unit cell is close 
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to unity, i.e., ria^ = « 1. The on-site en¬ 

ergy of conduction-electron orbitals fixes the energy zero: 
Ec = 0. 

Our main result is the phase diagram for the Ander¬ 
son model on the triangular lattice as obtained by site- 
dependent DMFT. This displayed in Fig. To cover 
the relevant parameter region, we have performed ~ 500 
independent DMFT calculations on the SuperMUC su¬ 
percomputer cluster of the LRZ Munich for different U 
and ^ in several massively parallel runs with step sizes 
AU = 0.5 and Afi — 0.05. We have considered the 
model, Eq. Q on a lattice with 25 x 25 unit cells and 
periodic boundary conditions to perform the fc-sum in 
Eq. 0 explicitly. This is sufficient to ensure that the 
results do not depend significantly on the system size as 
has been checked carefully. Self-consistent results are in¬ 
dicated as dots and symbols in Fig. in the U-n plane. 
About 200 iterations of the DMFT self-consistency cy¬ 
cle usually turned out to be sufficient for convergence. 
To allow for spontaneous breaking of the SU(2) spin- 
rotation symmetry, we explicitly treat the a =t and 
the a = 4 , channels as independent of each other within 
the CT-QMC solver. Furthermore, the DMFT cycle is 
started with a spin-asymmetric Hartree-Fock-type initial 
self-energy. It turns out that magnetic phases, if present, 
are easily found and stabilized in this way. Within the 
present study we focus on magnetic phases with collinear 
moments for simplicity even though non-collinear mag¬ 
netic phases may be expected in the case of the triangular 
lattice due to geometrical frustration. In fact, previous 



CO 


FIG. 2: (Color online) /- and conduction-electron densi¬ 
ties of states and respectively, for the non¬ 

interacting {U — 0) Anderson model on the triangular lattice. 
Energy units are fixed by the nearest-neighbor hopping t = 1. 
Further parameters: V = 1, Ec = £/ = 0. Centers of gravity 
are located at ca = 0. 


Hartree-Fock (HF) calculations at and off half-fillin^l^li^ 
suggest that a “classical” non-collinear 120° antiferro¬ 
magnetic phase is realized in a certain range of the phase 
diagram. We expect that, by enforcing collinearity, the 
120° phase is replaced by a collinear “titii” antiferro- 
magne tic p hase which has also been found within HF 
theory.E^ES 


A. Phase diagram 


Fig-i shows five different phases. At half-filling, the 
system is a non-magnetic Kondo insulator (KI) in the en¬ 
tire U range considered here. For fillings sightly off half¬ 
filling, the system stays non-magnetic but immediately 
becomes metallic. Above half-filling, this non-magnetic 
“Kondo singlet” (KS) phase is the only phase that has 
been found, at least up to n = 1.1 —1.2. Below half-filling 
and for a sufficiently strong interaction [/ > C4 « 2, there 
are two different magnetic phases, an antiferromagnetic 
phase (AFM) and a phase with partial Kondo screen¬ 
ing (PKS). The AFM phase is a collinear “t, tii” phase 
where the magnetic moments at two sites (say, A and 
B) in the unit cell are ferromagnetically aligned and of 
equal magnitude while the third moment is antiferro- 
magnetically oriented to the former two with a magni¬ 
tude such that the total magnetic moment in the unit 
cell is zero: rri'P + + mp = mg > 0 and 


Af) 


mP = —2mo < 0. Here, rria^ = n''A — and 


(/) 


Af) 


Af) 


at 


''at 


(c) _ (c) (c) 

m),’ = 

The PKS phase is characterized by one site (say A) 
with vanishing ordered magnetic moment, or almost van¬ 
ishing moment (see discussion below), while the mo¬ 
ment on the two remaining sites are of equal magnitude 
but antiferromagnetically aligned: mP + mP = mg = 
— (mP + mP) > 0. The total moment in a unit cell is 
again zero. The AFM and the PKS phases appear in a 
certain filling range UciiU) < nc 2 {U) which increases in 
width with increasing U and which is roughly centered 
around n « 0.9. The PKS phase appears at weaker U 
as compared to the AFM phase and separates the latter 
from the non-magnetic KS phase for n —>■ riciiU). At 
much lower fillings, there is also a ferromagnetic phase 
(FM) with a non-zero total magnetic moment per unit 
cell. This requires a significantly weaker critical interac¬ 
tion Uc « 1.25 as compared to AFM and PKS magnetic 
phases. The FM phase is realized in a rather narrow fill¬ 
ing range, roughly centered around n « 0.75 for weak U 
and n « 0.67 for U = i. 

We expect that the phase diagram obtained for inverse 
temperature /3 = 100 and shown in Fig. is close to 
the zero-temperature phase diagram. To estimate the 
remaining effects that are due to a finite /3, we have 
studied the parameter region close to the PKS phase for 
a somewhat higher temperature {/3 = 70). The results 
are shown in Fig. Comparing the phase diagrams for 
the different temperatures, there are no qualitative dif- 
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FIG. 3: (Color online) U-YS.-n phase diagram of the Anderson 
model on the triangular lattice as obtained by site-dependent 
dynamical mean-field theory. Each point corresponds to a 
converged DMFT calculation using CT-QMC (hybridization 
expansion, segment code) as a solver at /3 = 100. At half¬ 
filling n = 1 (solid line) the system is a Kondo insulator (KI) 
for all U > 0. Off half-filling, we find a metallic Kondo singlet 
state (KS, dots) as well as three different collinear magnetic 
phases: a partial Kondo-singlet phase (PKS, circles), an an¬ 
tiferromagnetic phase (AFM, squares) as well as a ferromag¬ 
netic phase (FM, triangles). 


ferences. Merely the extension of the AFM and the PKS 
phases in the U-n plane is somewhat reduced for /3 = 70, 
and the critical interaction increases a bit from [4 ~ 2 
(/3 = 100) to Ua « 2.5 (4 = 70). 
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FIG. 4: (Color online) U-vs.-n phase diagram as in Fig. [^but 
for /3 = 70 and in a smaller filling range including the AFM 
and PKS phases. 


B. Kondo insulator at half-filling 

We start the discussion with the KI phase at half¬ 
filling. The insulating nature of this phase is easily ver¬ 
ified by means of the charge susceptibility k = dn/dfi 
which is found to vanish at half-filling for any U > 0. 
For U = 0 and half-filling the system is actually a simple 
band insulator: The chemical potential is located in the 
hybridization band gap which opens for any K > 0, see 
Fig.H For the correlated system at U = 2.5, the charge 
gap Ac at half-filling can be read off from the /i range 
in which the charge susceptibility k vanishes; see lower 
panel of Fig. where is plotted as a function of the 
n. The gap persists for all C/ > 0 and decreases with in¬ 
creasing U as is obvious when comparing with the charge 
gap for U = 3.5, for example, which can be read off from 
the lower panel in Fig. (note the different scales for /i 
in the two figures). 

It is tempting to relate this decrease of the energy 
scale with increasing U to the decrease of the coupling 
constant J = 8V'^/U in the effective low-energy Kondo 
lattice that is formally obtained by the Schrieffer-Wolf 
transformatiorP in the local-moment regime of the An¬ 
derson lattice model. Local magnetic moments, required 
for magnetic long-range order, are formed on the / or¬ 
bitals due to a strongly repulsive Hubbard-C/. One must 
be aware, however, that even for U = A there are still 
substantial charge fluctuations on the / orbitals. This is 
indicated, for example, by a ^ 5% deviation of the av¬ 
erage / occupancy from unity at half-filling (see upper 
panels of Figs. [^and[^. Hence, the system is not fully 
in the local-moment limit. Nevertheless, we find an anti¬ 
ferromagnetic linear response of the conduction-electron 
magnetic moments when applying a homogeneous static 
magnetic field to the / electron spins. This indicates 
an antiferromagnetic (J > 0) coupling between the local 
/ and c spins consistent with the local-moment picture 
provided by an effective Kondo lattice. 

Deep in the local-moment regime for [/ —>■ oo at fixed 
V, the physics would be governed by a small energy scale, 
set by J, or even by Tk oc which makes calcula¬ 

tions at stronger U extremely difficult: In fact, we have 
not been able to stabilize a self-consistent solution of the 
DMFT equations at interaction strengths substantially 
stronger than {7 = 4. 

Interestingly, there is no magnetic phase found at half¬ 
filling n = 1. This is opposed to static mean-field (HF) 
theory for the same modeP^ which generates a rather 
complex phase diagram which comprises different mag¬ 
netic as well as insulating and metallic phases at half¬ 
filling. Using the Hartree-Fock approximation for the 
self-energy, 

= (7) 

we have reproduced the HF results of Ref. |T3 for V = 
1 as a check of our numerical implementation. As the 
DMFT correctly accounts for local fluctuations beyond 
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FIG. 5: (Color online) Upper panel: Difference between the 
occupancy of the / (c) orbitals and the average filling, —n 
{n^a^ — n), as function of n. Results for different sites in the 
unit cell: A (green), B (blue), C (red). Lower panel: Chemical 
potential/i as function of n. Calculations for f/ = 2.5. Further 
parameters: V = 1, P = 100. 


C. Metallic Kondo singlet phase 

For fillings slightly off half-filling, the system becomes 
immediately metallic and has a finite charge compress¬ 
ibility K > 0 (see lower panels of Figs. and |^. Actu¬ 
ally K turns out as non-zero for any hlling. Opposed to 
previous HF calculations,this implies that there is no 
instability towards phase separation. 

The local correlations between / and c spins are 
strongly antiferromagnetic as indicated by a correspond¬ 
ing antiferromagnetic linear response. Still, there is 
no magnetic order for hllings n < 1. We refer to 
this paramagnetic metallic state with local antiferromag¬ 
netic correlations as a heavy-fermion or Kondo-singlet 
state (KS) even if the local spin on the / orbitals, 

cannot be Seen as a rigid 
spin-S = 1/2 since the local / moment somewhat 

deviates from S{S -I-1) = 3/4. 

The respective top panels of Figs. and show the / 
orbital double occupancy relative to its non-interacting 
value, i.e., While the 

double occupancy is suppressed considerably for hllings 
close to half-hlling, it is still far from zero even at C7 = 3.5 
(Fig.§, for example, where Da ~ 0.45. At U = 2.5 (Fig. 

we hnd Da at a minimum for n « 0.65. 


the static mean-held theory, we conclude that those local 
fluctuations are sufficient to destroy any magnetic order 
at n = 1 (and in the U range considered here). 


0.85 


0.90 


0.95 
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FIG. 6: (Color online) The same as Fig. but for U = 3.5. 
Results for different sites in the unit cell: A (green), B (blue), 
C (red), as indicated. 


D. Antiferromagnetism 

Magnetic phases first appear at hllings centered around 
n « 0.92 for U = 2.5 and for U = 3.5 (Figs. and 
This is the hlling range where the / occupancy is at or 
very close to unity and where, despite substantial charge 
huctuations, the local-moment picture is most adequate. 
The magnetic coupling between the local moments must 
be provided by the a priori uncorrelated c orbitals, sim¬ 
ilar to the standard RKKY mechanisnP that can be de¬ 
rived perturbatively in the Kondo-lattice model. 

On the triangular lattice, however, magnetic order 
induced by indirect antiferromagnetic exchange is frus¬ 
trated. Except for a non-magnetic state, there are two 
obvious possible compromises to form a state with van¬ 
ishing total magnetic moment in the unit cell, namely a 
state with 120° orientations between pairs of magnetic 
moments as well as a collinear “t, ti i” phase. Apart 
from the PKS phase to be discussed below, the latter is 
the only plausible antiferromagnetic state if collinearity 
between the moments is enforced as is done here. 

For U = 3.5 and with decreasing filling n the sys¬ 
tem undergoes a phase transition to the AFM phase at 
n « 0.97. Fig. [^demonstrates that this phase transition 
is continuous with the staggered magnetization itiq (see 
definition above) as an order parameter that evolves from 
mo = 0 and increases with decreasing n in a continuous 
way. The magnetism is predominantly carried by the / 
moments with a maximum of Im^'^j « 0.6 while the c 
orbitals are by about one order of magnitude less polar- 
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FIG. 7: (Color online) Top panel: Filling dependence of the 
/-orbital double occupancy da'^ = (relative 

to its non-interacting value for U = 2.5. Results 

for different sites A (green), B (blue), C (red) in the unit 
cell as indicated. Middle and bottom panels: f and c ordered 
magnetic moments, and at the different sites in the 
unit cell for t/ = 2.5 as functions of n. 

ized (note the different scales in Fig. [^. Note that the 
site-dependent moments are oriented antiparallel to the 
respective / moments. 

Across the transition to the AFM phase there is hardly 
any change of the double occupancy i-e., 

the increase of Hq, seen in Fig. (top panel) is mainly 
due to the polarization of the / orbital only. For the 
“I” site in the “f, t: i” state this effect is a bit stronger 
as its magnetic moment has the higher absolute value. 
The fact that the double occupancy and thus the size 
of the local / moment is basically unaffected, favors a 
picture of magnetic ordering of preformed local moments 
and is consistent with an RKKY-like indirect exchange 
mechanism in the local-moment regime of the Anderson 
lattice. 


E. Partial Kondo screening 

With further decreasing n at U = 3.5 there is another 
second-order phase transition from the AFM state to a 
phase with partial Kondo screening (PKS) (see Fig. |^. 
For U = 2.5 the PKS phase directly evolves from the KS 
through a second-order transition (see Fig. [^. In both 
cases, a Kondo-singlet formed at one site in the unit cell, 
say A, coexists with a non-local pair of antiferromagnet- 
ically coupled moments at the B and C sites. The total 
ordered moment in a unit cell is zero. Eventually, for fill¬ 


FIG. 8: (Color online) The same as Fig. but for U — 3.5 
and for a smaller filling range. 


ings n < rici « 0.88 at 17 = 2.5 and for n < rid ~ 0.82 at 
U = 3.5, the system returns to a paramagnetic KS state 
in another continuous phase transition. 

As compared to the AFM and also, at even lower fill¬ 
ings, to the FM phase, the numerical stabilization of a 
self-consistent PKS solution is most difficult, i.e., a large 
number of iterations (up to 200) is required. This also 
reflects itself in the remaining (unphysical) noise on the 
PKS data seen in Figs. [^and[^ As a technical remark, 
let us mention that each DMFT run is completely in¬ 
dependent from the preceding one and starts from the 
same initial guess for the self-energy which is taken as fre¬ 
quency independent and homogeneously spin-polarized. 
Due to this independency, the self-consistent values for 
the magnetic moments typically do not always form con¬ 
tinuous functions of /i, because arbitrary permutations 
of the A, B, C sites in a unit cell and also a global sign 
change a ^ —a yield physically equivalent solutions of 
the DMFT equations. We have employed those symme¬ 
try operations a posteriori in scans with extremely small 
steps in the chemical potential (A/r = 0.007) to generate 
functions as continuous as possible by means of least- 
square fits minimizing the parametric distance between 
pairs of consecutive self-consistent solutions. 

The Kondo effect requires a locally antiferromagnetic 
effective coupling between the local / and c spins. This is 
clearly present: As mentioned above, the linear response 
of the c moments to a static magnetic field applied to the 
/ electron spins is found as antiferromagnetic in param¬ 
agnetic phase close to the AFM and PKS phase. Fur¬ 
thermore, within the symmetry-broken PKS phase, the 
ordered moments rria'^ and rria'^ are antiferromagneti- 
cally aligned on the B and C sites. On the other hand. 
































the robustness of the PKS phase, i.e., its extension in 
the U — n plane, and also the presence of strong charge 
fluctuations, see the sizable double occupancy in Figs. 
and|^ suggest that the physics is non-universal and by 
no means ruled by a single Kondo scale Tk- 

It is interesting to note that our data unambiguously 
show that there is no “perfect” partial Kondo screening. 
Namely, a slight polarization < 0 and > 0 
of the local / and c spins on the A, i.e., on the Kondo 
site is clearly visible in Figs. and The proximity to 
the pair of RKKY-like antiferromagnetically coupled mo¬ 
ments, which explicitly breaks time-reversal symmetry, 
implies that there are admixtures of states with non-zero 
spin quantum number to the Kondo “singlet”. Assum¬ 
ing this admixture to be given by a single spin-triplet 
state for simplicity, the antiferromagnetic environment 
explains a coupling to the M = 0 component of the 
triplet. A finite polarization of the Kondo singlet, how¬ 
ever, rather requires a coupling to the M = ±1 compo¬ 
nents and thus implies the additional breaking of the Z2 
symmetry of the antiferromagnetic state. This sponta¬ 
neous symmetry breaking in the PKS phase is also visible 
in the magnitudes of the B- and C-site moments, namely 
Im^^l > and is present in the AFM state 

anyway. 

Accompanying the ordering of the spin degrees of free¬ 
dom, there is a also a (weak) charge ordering in the AFM 
and the PKS phase (see upper panels of Figs. [5 and[^: 
There are two interesting observations: First, the devi¬ 
ation of the charge density from the average density is 
much stronger on the c orbitals as must be expected in 
the local-moment regime where charge fluctuations on 
the correlated / orbitals are very effectively suppressed. 
This effect is stronger for U = 3.5 and compared to 
U = 2.5. Second, within the PKS phase, there is a 
charge transfer from the “Kondo site” (A) to the “mag¬ 
netic sites” (B, C): c- 

the Kondo effect, the local conduction-electron density of 
states at the A site will develop a dip, and spectral weight 
must be shifted above or below the Fermi energy. In the 
absence of particle-hole symmetry, this shift is asymmet¬ 
ric and changes the occupancy. The sign and the size 
of the resulting charge transfer, however, depend on the 
details of the band structure. Charge disproportiona¬ 
tion was also found within the PKS (“partial disor der”) 
state that is obtained by means of the HF approach. h^h^ l 
Opposed to our DMFT results, the charge transfer seen 
in the HF studies is much larger for the / as compared 
to the c orbitals. This must be seen as an artifact of 
the static mean-field approach which cannot account for 
local-moment formation. 

As a function of U, the PKS phase is located between 
the KS and the AFM phase in the phase diagram. This 
can be understood by referring to the famous Doniach 
diagram®! In the KS phase at weaker U (stronger J) 
the Kondo effect dominates while for strong U (weak J) 
the RKKY interaction is dominant and results in mag¬ 
netic order. The PKS state can be seen as a possible way 


to avoid geometrical frustration in the antiferromagneti¬ 
cally ordered state which is preferred if the formation of 
a Kondo singlet is less expensive than breaking up two 
frustrated magnetic bonds and forming a non-frustrated 
third one. As a compromise between indirect exchange, 
frustration and the Kondo effect, it appears between the 
KS and the AFM phase. 


F. Ferromagnetism 

At lower fillings around n = 0.7, depending slightly on 
U, the system develops homogeneous ferromagnetic order 
(see Fig. [^. As can bee seen from Fig. [^for U = 2.5, the 
transition to this state is continuous at the lower as well 
as at the upper critical density. The ferromagnetic state 
is metallic with a finite compressibility (see Fig. and 
partially polarized with a maximum ordered / moment of 
« 0.52 at n « 0.71. The moment on the conduction- 
electron orbitals {'rrS‘^'> « 0.02) is more than an order of 
magnitude smaller and ferromagnetically aligned to the 
moment on the / orbitals. 

Generally, there are several mechanisms that may 
cause metallic ferromagnetism®!!^ The main idea of the 
RKKY concept!^ consists in a magnetic coupling of well- 
formed local / moments in an effective Kondo-lattice 
modep! which is mediated by the conduction electrons 
and features ferromagnetic order if the effective RKKY 
coupling Jrkky(<7) = = 0 ) is peaked at 

q = 0. While the RKKY theory is a perturbative ap¬ 
proach (J —>■ 0), the double-exchange mechanisnPSEIl 
applies to the strong-J regime of a Kondo lattice and 
predicts that a ferromagnetic ordering of the / moments 
minimizes the kinetic energy of the conduction electrons. 

It is questionable, however, if those concepts apply here 
as there are strong charge fluctuations preventing the for¬ 
mation of well-defined / moments in our case. This is ob¬ 
vious from the sizable deviation of the / occupancy from 
unity (see Fig.« 0.8 in the relevant filling range). 
Another clear indication that the system is no longer in 
a local-moment regime is the ferro- rather than antiferro¬ 
magnetic coupling between / and c moments (see middle 
and lower panel of Fig. [^around n = 0.7). This is in¬ 
compatible with an effective low-energy Kondo model. 

It is interesting to note that this implies a filling- 
dependent crossover from the local-moment regime with 
a locally antiferromagnetic coupli^ between / and c mo¬ 
ments (see and m.a'^ in Fig. [^in the PKS and AFM 
phases) to a mixed-valence regime. This can also be ver¬ 
ified easily by studying the linear response in the param¬ 
agnetic phase separating the FM and the PKS phase in 
Fig.i By applying a weak magnetic field to the / mo¬ 
ments, one finds that the local coupling between / and c 
moments changes from antiferro- to ferromagnetic with 
decreasing filling. 

At {7 = 0 the static off-diagonal f-c magnetic spin 
susceptibility can be computed easily in the entire filling 
range. Except for low fillings around and below n « 0.25, 
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corresponding to the van Hove singularity of the density 
of states close to the lower band edge (see Fig. [^, the 
local response is found as antiferromagnetic for n < 1. 
Above half-filling, the response turns to ferromagnetic 
and is at a maximum for n « 1.2 corresponding to the 
van Hove singularity at w « 0.4 (see Fig.|^. We conclude 
that the ferromagnetic phase cannot be understood as 
an instability of the Fermi sea in the weak-t7 regime. 
Just the opposite, the paramagnetic state from which 
the ferromagnetic phase evolves should be considered as 
strongly correlated. Already for U = 2.5, the double 
occupancy is strongly suppressed and Da is in fact at a 
minimum for n « 0.65 (see Fig. [^. 

The importance of a strong asymmetry of the density 
of states for metallic ferromagnetism at strong and in¬ 
termediate interaction strengths has been empha sized by 
DMFT studies of the single-band Hubbard modeL^^^^^HE] 
The key idea is that in a situation where double occu¬ 
pancies are effectively suppressed, the system does not 
gain much interaction energy from ferromagnetic order¬ 
ing. Therefore, the appearance of ferromagnetism must 
be understood by referring to the (complicated) kinetic 
energy of the correlated paramagnetic state from which 
it derives. Within DMFT this suggests that the shape 
of the non-interacting density of states becomes impor¬ 
tant. In fact, studying the impact of a shape-controlling 
parameter ,22 ferromagnetism was demonstrated to be fa¬ 
vored in cases with a highly asymmetric density of states, 
in a parameter range where the density of states is high^ 
and at strong to intermediate interaction strengths. 

We propose that a similar line of reasoning applies to 
the periodic Anderson model in the considered parame¬ 
ter region: Even at U = 2.5 and all the more for stronger 
U, double occupancies are strongly suppressed, and the 
gain in kinetic energy obtained by ferromagnetic order¬ 
ing is dictated by a strongly asymmetric partial / den¬ 
sity of states. The filling range where ferromagnetism is 
likely to occur, is then indicated by a corresponding high 
density of states. Note that n = 0.7 corresponds to a 
non-interacting chemical potential of /i « —0.62 which is 
already close to the van Hove singularity (at uj = —0.18). 
Substantially higher fillings would be even more favor¬ 
able for ferromagnetism, but here the crossover to the 
local-moment regime and the developing antiferromag¬ 
netic correlations overwrite the ferromagnetic tendencies. 

This picture also explains why the FM phase shifts 
to lower fillings with increasing U in Fig. Stronger 
interactions favor ferromagnetism and extend the FM 
phase to a larger filling range as is again well known 
from the single-band case.l^ This explains the decrease 
of the lower critical filling for the FM phase with increas¬ 
ing U. At the same time, however, an increasing U favors 
local-moment formation, and therefore the KS phase with 
antiferromagnetic correlations extends at the cost of the 
mixed-valence regime. This explains the decrease of the 
upper critical filling with increasing U. 

Previous worlPsHUl on ferromagnetism in the periodic 
Anderson model has been done using different theoreti¬ 


cal approaches and in largely different parameter regimes. 
Nevertheless, ferromagnetic order away from half-filling 
appears as a robust result. As basically all studies have 
exclusively been performed for bipartite lattices, a direct 
comparison with our results is not possible. There are, 
however, close similarities with the results of a DMFT 
study by Meyer and NoltinJ^ which, for a Bethe lat¬ 
tice with infinite connectivity, demonstrates that ferro¬ 
magnetism appears in the mixed-valence regime for a fi¬ 
nite filling range. This study also points out a crossover 
from antiferro- to ferromagnetic coupling between / and 
c magnetic moments with decreasing filling, consistent 
with our findings, and suggests a mechanism based on an 
effective single-band model with strongly correlated and 
itinerant electrons - an idea that was formalized later on 
by Batista et al.^^ 


IV. CONCLUSIONS 

Our site-dependent DMFT study of the magnetic 
phase diagram of the periodic Anderson model on the 
triangular lattice has uncovered a surprisingly complex 
phenomenology which could be traced back to a compe¬ 
tition between several physical mechanisms at work. In 
particular, the phase diagram is governed by: 

(i) the formation of local magnetic moments on the / 
orbitals. Due to the non-bipartite structure of the tri¬ 
angular lattice, half-filling of the / orbitals is found for 
total fillings below half-filling, around n « 0.9, weakly 
depending on U. Here, the low-energy physics is well cap¬ 
tured by an effective Hondo lattice although there are siz¬ 
able / charge fluctuations for the weak- to intermediate¬ 
coupling regime considered here (D < 4). At somewhat 
lower fillings, there are still well-developed local / mo¬ 
ments, but the charge fluctuations increase since the / 
electrons become itinerant. 

(ii) Mixed-valence physics with strong charge fluctu¬ 
ations on the / orbitals, even at stronger U, replaces 
the local-moment regime for lower fillings (roughly be¬ 
low n « 0.75, depending on U). The filling-dependent 
crossover from the local-moment to the mixed-valence 
regime is accompanied by a reversal of the effective local 
exchange between the local / and c spins from antiferro¬ 
magnetic to ferromagnetic. 

(iii) An RKKY-like indirect magnetic exchange be¬ 
tween the / magnetic moments induces antiferromag¬ 
netic order for sufficiently strong U within the local- 
moment regime. As we have enforced spin structures 
to be collinear, this results in an “f, f,AFM phase on 
the triangular lattice which possibly mimics “classical” 
120° AFM order. 

(iv) The Kondo effect competes with the indirect ex¬ 
change in the spirit of the Doniach diagram. At low tem¬ 
peratures, besides magnetic ordering, the large entropy 
carried by the local-moment system can be removed by 
screening the / moment in a Kondo singlet with the 
conduction-electron spin degrees of freedom. With de- 
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creasing U, and prior to charge fluctuations becoming 
dominant, this Kondo-singlet (KS) phase replaces the an¬ 
tiferromagnetic order. Kondo physics is also dominating 
for lower fillings around n ^ 0.8, depending slightly on 
[7, as well as for fillings very close to and at half-filling. 
The hybridization band gap in the non-interacting den¬ 
sity of states results in a band insulator at half-filling for 
[7 = 0 and develops into a correlated Kondo insulator 
with increasing U. 

(v) Geometrical frustration affects the competition be¬ 
tween Kondo screening and RKKY coupling. At the bor¬ 
der between the AFM and KS phase, it becomes favor¬ 
able to avoid frustration by partial Kondo screening of 
one / moment per unit cell. This allows the remnant mo¬ 
ments to form an unfrustrated RKKY-coupled collinear 
antiferromagnet. The PKS phase is metallic, and it sup¬ 
ports a (weak) charge-density-wave ordering in addition, 
mainly on the c orbitals. Although it results from a 
compromise between Kondo screening, RKKY coupling 
and frustration, the PKS state has turned out to be sur¬ 
prisingly robust. It appears in an extended parameter 
range and does not need any anisotropic terms in the 
Hamiltonian.lS! Due to proximity to the time-reversal- 
symmetry-breaking RKKY-coupled remnant moments, 
the partial Kondo screening is imperfect resulting in a 
tiny magnetic moment on the / and, antiferromagneti- 
cally aligned, on the c orbital at the “Kondo site”. 

(vi) Strong correlations among itinerant electrons give 
rise to a metallic and partially polarized ferromagnetic 
phase in the mixed-valence regime. In this case, the 
non-bipartite lattice structure favors magnetic order as 
it produces a highly asymmetric non-interacting density 
of states which is known to crucially affect the kinetic- 
energy balance favoring ferromagnetism in a range of fill¬ 
ings with high density of states at the Fermi level and 
where antiferromagnetic correlations are subdominant. 
As a non-perturbative effect, this itinerant-electron fer¬ 
romagnetism lacks a clear (simple) mechanism - even in 
a single-band Hubbard model. 

Two main results of our study might be relevant for 
the understanding of PKS in real materials, such as 
CePdAl,E21 UNi 4 l^^ or even artificial geometries of mag¬ 
netic atoms on metallic surfaces,!^ and for corresponding 
electronic-structure models: (i) The PKS state appears 
at non-integer fillings. One might thus speculate that 
the gain in kinetic energy is essential to stabilize the 
state and that spin-only models may be questionable, 
(ii) The PKS state exclusively shows up at the border 
between the paramagnetic heavy-fermion and the mag¬ 


netically ordered phase. This could be tested experimen¬ 
tally by steering the system through this border, either 
by controlling the temperature or by means of chemical 
substitution.ES 

There are several lines along which our study could be 
continued in the future: First, non-collinear phases are in 
principle accessible by an SU(2)-symmetric formulation 
of the site-dependent DMFT. This may lead to a cer¬ 
tain refinement of the magnetic phase diagram, with non- 
collinear (or even incommensurate) AFM phases partially 
replacing the “f, ti i” phase but we do not expect a fur¬ 
ther qualitative change as the relevant energy scale is still 
set by the effective RKKY-exchange coupling constant. 

Second, it would be interesting to make contact with 
the corresponding phase diagram of the Kondo model 
on the triangular lattice, either by applying DMFT to 
the Kondo model directlji^ or by using a solver which 
allows to resolve the Kondo scale, such as the numerical 
renormalization group.^ 

Finally, one may address the effect of non-local corre¬ 
lations beyond the single-site DMFT. The Kondo effect 
results from correlations between a single correlated / 
orbital and the conduction-band system and is, there¬ 
fore, captured correctly by a dynamical mean-field the¬ 
ory which treats those correlations exactly. DMFT also 
provides an accurate description of the non-local indirect 
exchange but the feedback of non-local magnetic correla¬ 
tions on the self-energy is missing. Those missing fluctu¬ 
ations must result in mean-field artifacts. Typically, the 
(site-dependent) DMFT approach is therefore, to some 
extent, biased towards the formation of magnetic order 
and tends to favor a symmetry-broken state | t) I i) at the 
expense of a non-local singlet (| t) I i) — I i) I t))/ 

One might speculate that, compared to the PKS state, 
the AFM phase is overestimated and that both, the PKS 
and the AFM phases are overestimated as compared to 
the KS state. 
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